Spontaneous spin current near the interface between unconventional superconductors 

and ferromagnets 



^1- 
O 

o 

> 
O 



o 
o 

c3 



I 

C 

O 

o 



> 

(N 

m 

O 

-i— > 
c3 



I 

c 

o 
o 



Kazuhiro Kuboki and Hidenori Takahashi 
Department of Physics, Kobe University, Kobe 657-8501, Japan 
(February 2, 2008) 

We study theoretically the proximity effect between ferromagnets (F) and superconductors (S) 
with broken time-reversal symmetry (T). A chiral (p x ±ip H )-wave, and a d x 2 _ y 2-wme superconduc- 
tor, the latter of which can form T-breaking surface state, i.e., (d x 2_ y 2 ±zs)-state, are considered for 
the S side. The spatial variations of the superconducting order parameters and the magnetization 
are determined by solving the Bogoliubov de Gennes equation. In the case of a chiral (p x ± ip y )- 
wave superconductor, a spontaneous spin current flows along the interface, but not in the case of a 
da;2_ H 2-wave superconductor. For F/S(p x ± ip y )/F trilayer system, total spin current can be finite 
while total charge current vanishes, if the magnetization of two F layers are antiparallel. 



PACS numbers: 74.45. +c, 74.50.+r, 74.78. Fk 



I. INTRODUCTION 

Recently the proximity effect of unconventional super- 
conductors has been a subject of intensive study [1-8]. 
This is because the interface properties of these supercon- 
ductors can be quite different from those of conventional 
(s-wavc) ones due to the nontrivial angular structure of 
pair wave functions, so that their study is of particular 
interest. The proximity effect between unconventional 
superconductors and magnetic materials is an important 
problem, since there are competition and the possible co- 
existence of two kinds of orders. The unconventional su- 
perconductivity is usually induced by the magnetic inter- 
action, and the existence of nodes in the gap is favorable 
for the coexistence compared with the case of conven- 
tional superconductors. If the magnetism and supercon- 
ductivity coexist near the interface, the electronic prop- 
erties of the surface state become quite unusual. In par- 
ticular the spin-triplet superconducting order parameter 
(SCOP) can be induced near the surface of spin-singlet 
superconductors [5]. (Similar effects for conventional s- 
wave superconductors and ferromagnets have been stud- 
ied in [9-11].) 

In this paper we study the proximity effect between 
ferromagnets and unconventional superconductors with 
dj;2_ y 2-wave and (p x ± ipj,)-wave symmetries. The for- 
mer state is realized in high-T c cuprates [12], while the 
latter is a candidate of the superconducting (SC) state 
of Sr2Ru04 [13]. It is known that near the [110] sur- 
face of the d a .2_ !/ 2-wave superconductor the system may 
break time-reversal symmetry (T) by introducing second 
component of SCOP [14-18,7,8], and that a spontaneous 
current and fractional vortices may occur at the surface 
if T is broken [19-21,14-18,7,8]. The {p x ± ip y )-wave SC 
states also break T and the spontaneous current arises 
at the edge of the system [22]. When a ferromagnet is 
attached to these unconventional superconductors, the 
magnetization may be induced in the latter due to the 



proximity effect [5]. Then it may be expected that a 
spontaneous spin current can appear along the interface 
between a T-breaking superconductor and a ferromagnet, 
because there is an imbalance of the densities of spin-up 
and spin-down electrons. We will show that it is actually 
possible for the superconductor with (p x ±ip y )-wave sym- 
metry, but not in the case of d x 2_ y 2-w&ve symmetry [23]. 
(Spin currents in the case of interfaces between s-wave 
superconductors and ferromagnets have been found in 
[24], though their treatment of the magnetization is not 
fully self-consistent.) 



II. MODEL AND BDG EQUATIONS 

The system we consider is a two-dimensional fer- 
romagnet(F)/superconductor(S) bilayer and F/S/F tri- 
layer systems, and we treat [100] and [110] interfaces as 
schematically depicted in Fig.l. The directions perpen- 
dicular (parallel) to the interfaces are denoted as x (y) 
and x' (y'), for [100] and [110] interfaces, respectively. 
The axes x' and y' are 45° rotated from the crystal axes 
directions x and y. We assume that the system is uniform 
along the direction parallel to the interface. In order to 
describe the magnetism and superconductivity, the tight- 
binding model on a square lattice with on-site repulsive 
and nearest-neighbor attractive interactions is treated 
within the mean-field (MF) approximation [5,6], and we 
consider only the case of zero temperature (T = 0). The 
Hamiltonians for the two layers are given by 



h l = -t L ( c U« 



■3," 



n iT n a 



<IJ>CT 



(i) 



<i,3> 



where denotes the nearest-neighbor bonds, 1%^ = 

cl a Ci a , and Cj CT is the annihilation operator of electrons 
at a site i with spin a. Parameters tj,, Ul and Vl 
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represent the transfer integral, the on-site interaction 
and the nearest-neighbor interaction, respectively, for the 
L(= F,S) side. The transmission of electrons at the in- 
terface is described by the following tight-binding Hamil- 
tonian 



H T = -t T E ( C l,a C m,<r + h.C.) 



(2) 
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where I (m) denotes the surface sites of F (S) layer, 
and then the total Hamiltonian of the system is H = 
Hp + Hs + Ht — /i c t<j c i<? with \x being the chemical 
potential. The interaction terms are decoupled within 
the MF approximation: 

Un^n^-f U{n^)n li + U{n l[ )n l] -U{n^){n l < l ), 
Vn nnjl VAycJ^ + VA* ]Cl1 c n - V\A^ [d) 

with Ajj = {ctfCji). Then Ay and the magnetization 
m, = (n,f — nii)/2 (perpendicular to the plane) are 
the OP's to be determined self-consistcntly. The SCOP 
with each symmetry can be formed by combining Ay's: 
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A {a > - + A (s > + A (5) 0/4 
(extended s-wave ), A px(y) (i) = (Ag)_ s(fl) - Agl £( . ) )/2 

( S ) = ( A . . _i_ A .A /O o„^l A( T ) = 



( Px(y) -wave), where A^ ; = (Ay + A^)/2 and A^ > = 
(Ay — Ajj)/2 are the spin-singlet and the spin-triplet 
pairing OP's, respectively. 

We impose the open (periodic) boundary condition for 
the x and x' (y and y') directions for [100] and [110] 
interfaces, respectively, and carry out the Fourier trans- 
formation along the y- and y '-directions. For the [100] 
case we define A^ = Ai^± x and A y± = Ai_i± y as the 
SCOP's independent of y. Similarly for the [110] case, 
A? = A M+Q ^ (a, = + or -) are SCOP's inde- 
pendent of y' . (\x\ = \y\ — a, \x'\ = \y'\ = a' with a and 
a' defined in Fig.l.) Then the mean-field Hamiltonian is 
written as (hereafter i denotes Xi and x i for a [100] and 
a [110] interface, respectively) 

Wmfa = E E E *!WM*0*;(*0 (4) 

k i j 

with ty\(k) — (c\^(k),Cii(— k)), and k is the wave num- 
ber for the direction parallel to the interface. The matrix 
hij(k) is given as 
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where 
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for the [100] interface, and 

iija{k) = {-/j, + U(n i - a ))S ij 

_ 2t ij cos(— )(5. +6.. ./) 

F«(fc) = Vtfp. ._ ^(At+e ifea '/2 + A +- e -ifca'/ 2) 



(7) 



+ 5. . + ^(Ar+ e ^72 + A - 



3 — ifca'/2\ 



for the [110] interface. Here iy — tp(s)> Vij = Vf(S) if 
both i and j are on the F (S) side, while ty = tx, V.j = 
if i and j correspond to the interface sites. 

We diagonalize the mean-field Hamiltonian by solving 
the following Bogoliubov de Gennes (BdG) equation [25]: 



EM fc )(";; (fc) 



Uin {k) 



where E n (k) and («,„(!;), Vi n (k)) are the energy eigen- 
value and the corresponding eigenfunction, respec- 
tively, for each k. The unitary transformation using 
(ui n (k),Vi n (k)) diagonalizes the matrix TCmfa, and con- 
versely the OP's Aij and rrii can be written in terms of 
E n (k) and (u in (k),v in (k)): 

y k,n 

A * ± = ^Y, u ^ k Kn(k)[l-f(E n (k))}e Tika 

y k,n 

Af = ^7 E UiAk)v* +a% ,Jk)[l f{E n {k)))e-^ ka '^ 

y k,n 

^T)= ]J-El«i,n(A)| 2 /(^n(fc)) 

(9) 



where /(i?) is the Fermi distribution function. (For the 
[110] case, N y in the last two expressions should be re- 
placed by N y t.) These constitute the self-consistency 
equations which will be solved numerically in the follow- 
ing. 

The procedure of the self-consistent numerical calcula- 
tion is the following. We substitute an initial set of OP's 
in the matrix elements of hij(k), and solve BdG equation 
(Eq.(8)) to get eigenvalues and eigenfunctions. Then we 
recalculate the OP's, and the iteration is performed until 
the values of OP's are converged. We have used various 
sets of initial OP's for the same values of the parameters 
of the model. If several different solutions are obtained, 
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we adopt the solution with the lowest energy as the true 
one. 

In the uniform case with U > and V = (i.e., repul- 
sive Hubbard model), the ground state phase diagram 
within the MF approximation was examined by Hirsch 
[26]. We use this result to choose the parameters to re- 
alize the ferromagnetic state in the F layer. The phase 
diagram in the case of U = and V < was also stud- 
ied, and the d x 2_ y 2-, extended s- and (p x ± ip y )-w&ve SC 
state appears depending on the band filling [27]. Various 
SC states can occur in the model with a single type of 
(nearest-neighbor attractive) interaction because of the 
change of the shape of the Fermi surface [28] . Near half- 
filling (fi ~ 0) the d x 2_ y 2-wa,ve SC state is stabilized, 
while the extended s-wave SC state is favored near the 
band edge (/i ~ ±4i). In the region between d- and s- 
wave states spin-triplet (p x ±ip y )-w&ve SC state appears. 

The properties of the interface states do not depend 
on the values of the parameters in a qualitative way, un- 
less the symmetries of the superconducting (d x 2_ y 2 or 
(p x ± ip y )-wave) and the magnetic (ferromagnetic or an- 
tiferromagnetic) states are changed. Therefore we will 
choose typical values of the parameters in order to real- 
ize the states in question. We will use finite values of Up 
to investigate the effect of electron correlations in the 
S side, while the attractive interaction in the F side is 
assumed to be absent, i.e., Vp — 0. Hereafter we take 
tp = ts(= i) = 1 as the unit of energy, and the tunneling 
matrix element tx is varied to see how the extent of the 
proximity effect is changed. 

III. INTERFACE STATES AND SPONTANEOUS 
SPIN CURRENTS FOR CHIRAL P-WAVE 
SUPERCONDUCTORS 

In this section we study the proximity effect between 
ferromagnets and superconductors with (p x ± ip y )-wave 
symmetry. Here the parameters are chosen to realize 
the ferromagnetic and (p x ± ip y )-w&ve SC states in F 
and S layers, respectively, i.e., Up = 12, Vs = —2.5 and 
l-i = —1.6. The system size we used is N x = N y = 120, 
which is large enough to study the interface states, be- 
cause the coherence lengths of SC and ferromagnetic 
states are of the order of 10a for the parameters used 
here. The spatial variations of the SCOP's and the mag- 
netization to near the [100] interface arc shown in Fig. 2 
where the tunneling matrix element is taken to be tx = 1 , 
and Us = 1.5. It is seen that A px and A py are suppressed 
near the interface, and the spin-singlet components 
and A s are induced. Near a surface of a chiral supercon- 
ductor faced to vacuum, the spin-singlet components are 
not induced, though the p-wave SCOP's are suppressed 
also in this case [22]. This difference is due to the pres- 
ence of to in the S side, because the SC state cannot be 
formed with (p x ± ip y )-channel only, if (nj) ^ (ri|). The 



proximity effects will be reduced as the tunneling ma- 
trix element becomes smaller. This is actually the case 
as seen in Fig. 3 where tx = 0.1 is used. It is seen that 
the penetration of to into the S side is suppressed, and 
d- and s-wave SCOP's are much smaller compared with 
Fig. 2. The change of the surface states, as tx is varied, is 
rather simple in the case of (p x ±ip y )-wa,ve superconduc- 
tors. Namely, the state monotonously approaches that 
for tx — as tx is reduced. (On the contrary, the surface 
state of a d-wave superconductor shows a phase transi- 
tion when tx is varied, as we will see in the next section.) 
The penetration depth of to depends on the electron cor- 
relation in the S side. The spatial variations of OP's for 
Us = 5 and tx = 1 are shown in Fig. 4. For this value 
of Us the system gets closer to the magnetic instability 
(though this Us is not large enough to stabilize the mag- 
netic order), and the correlation length of to becomes 
larger. Then the magnetization to penetrates far inside 
the S side compared with Fig. 2, and the d- and s-wave 
SCOP's decay rather slowly. 

In the T-violating SC state a spontaneous current flows 
along the surface in the region where the SCOP's are 
spatially varying [19-21]. Note that the condition for the 
minimum free energy requires the current normal to the 
interface to vanish. This condition in turn leads to a fi- 
nite current along the interface if T is broken. (Namely 
the present spontaneous current is the equilibrium cur- 
rent.) The expressions for the currents carried by spin-up 
and spin-down electrons are given by 

4 (i) = -^^E sin ka\ Ul Ak)\ 2 f(E n (k)) 

y k,n 

4» = ^^E sinfca K«(fc)l 2 [i-/(^(fc))]- 

y k,n 

(10) 

Since we do not consider the effect of the vector poten- 
tial A in this paper, terms linear in A in the expres- 
sions of J y are not taken into account. Near the sur- 
face faced to vacuum only the charge current (defined as 
-^charge = J T + J^) can flow [22]. When the surface is at- 
tached to a ferromagnet, to is finite in the region where 
the chiral SCOP's have spatial variations. Thus is not 
equal to J- 1 , and the spontaneous spin current (defined as 
J sv \n = J T — J^) appears as seen in Fig. 5 (a). For small 
value of tx(= 0.1), only small spin current flows while 
the charge current is large (Fig5(b)). This is because the 
penetration of to into the S side is suppressed (Fig. 3), and 
then the difference between Jj and J| is reduced. When 
the (p x ± zp 3/ )-wave superconductor is faced to vacuum, 
only a charge current arises and no spin current shows up 
[22]. The ir-dependence of spin currents presented here 
is consistent with this fact. When the effect of electron 
correlations is important in the S side, J c harge is reduced 
but J sp in is not (Fig. 5(c)). This is because a large Us 
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suppresses A p 's leading to a reduced J c hargc, while it en- 
hances m and J sp m- The spin current also becomes larger 
as the value of Uf is increased (Fig. 5(d) with Up = 16), 
since m in the F side is enhanced leading to larger value 
of m in the S side. To summarize the discussion of the 
spin current in bilayer systems, J sp i n can be large if the 
magnetization induced in the S side is large. 

We have also investigated the case of [110] interface 
between a (p x ± ip y )-wave superconductor and a ferro- 
magnet, and the results are qualitatively the same as 
those for the [100] interface. 

Now we examine the F/S/F trilayer system where S 
is a chiral (p x ± ip y )-wave superconductor. The spa- 
tial variations of the magnetization m and the p-wave 
SCOP's are shown in Fig. 6. (Small ReA^ and ReA s are 
also induced near the interface as in F/S bilayer system, 
though they are not shown here.) In the F/S bilayer sys- 
tem, we have seen that spontaneous spin currents as well 
as charge currents flow along the interface. The direc- 
tion of the charge current, J c hargo, is determined by the 
chirality (i.e., (p x + ip y ) or (p x — ip y ) ), and the currents 
at the opposite edges of the superconductor should flow 
in opposite directions. The direction of the spin current, 
^spin, is the same as (opposite to) that at the other edge 
of the superconductor if the magnetizations of two fer- 
romagnets are antiparallel (parallel). Then if the chiral 
superconductor is sandwiched by ferromagnets with an- 
tiparallel magnetizations, the total charge current (inte- 
grated over x) should vanish, while the total spin current 
will be finite. On the contrary if the magnetization of 
two F layers are parallel, both the total charge current 
and the total spin current should vanish. This is actually 
the case as shown in Fig. 7, where J2i <^chargc(«) vanishes 
for both parallel and antiparallel configurations of m and 
J2i Jspin{i) is finite only for an antiparallel configuration 
of m. If we assume t = leV, J2i ^s P in(«) ~ 10/iA. 



magnetic field [30,31].) In the (<i±is)-wave surface state, 
the spontaneous charge current flows along the y' direc- 
tion. If a ferromagnet is attached to a [110] surface of a 
dj;2_ y 2-wave superconductor, the coexistence of the mag- 
netization m in the T-violating state may be expected, 
leading to a spontaneous spin current. 

In Fig. 8 the spatial variations of m and the SCOP's 
near the [110] interface are shown for tx = 0.1. Since 
a d x 2_ y 2-wa,ve SC state has nodes in contrast to a fully 
gapped (p x ± ip y )-wave SC state, the proximity effect 
in the former case is much larger than that in the lat- 
ter. (See Fig. 3 where the same tx = 0.1 is employed in 
the case of (p x ± ip y )-wave superconductor.) When tx is 
not so small (tx > 0.02 for U's and V's used here), the 
magnetization penetrates into the S layer (Fig. 9). When 
m^0, the density of spin-up and spin-down electrons are 
not equal, and the Cooper pairs cannot be formed in sin- 
glet channels only, and thus the p-wave components have 
to arise [5]. We find that the complex component (is), 
which could be induced near a surface faced to vacuum, is 
destroyed, and the resulting state has a (d x 2_ y 2+p x +p y )- 
symmetry. When the value of tx is reduced, the symme- 
try of the SC state changes. For tx ^5 0.02, (d±is)-wave 
state appears, where the magnetization m (and thus p- 
wave SCOP's) is absent in the S layer even though tx 7^ 0. 
(Fig. 10) 

In the (d ± zs)-state a spontaneous charge current can 
flow along the interface, but there is no spin current be- 
cause m = 0. In the (d + p x + p y )-state both charge 
and spin current vanish, because all SCOP's are real. 
(Fig.ll) For the [110] interface the expression for the 
current J y/ (i) is slightly different from that for [100] 
case. We checked that the current from the site i to 
i + x' /2 + y'/2, denoted as J„ (i), is the same as that 
from the site i + x'/2 + y'/2 to i + y', denoted as J£(i). 
(Otherwise J x ,(i) ^ 0.) Then J y ,(i) is given by the y' 
component of J° (i) (or, equivalently J%(i))\ 



IV. INTERFACE STATES FOR D-WAVE 
SUPERCONDUCTORS 

In this section we examine the interface states for fer- 
romagnets and superconductors with d x 2_ y 2-w&\e sym- 
metry. In order to realize the ferromagnetic and the 
d x 2_ y 2-w&ve SC states in F and S layers, respectively, 
the parameters are chosen to be Uf — 10, Us = 1.5, 
Vs = —1.5 and fi = —0.5. We consider only the [110] 
interface here, since a T-breaking surface state may be 
formed in this case, but not in the case of a [100] sur- 
face [29]. The T-violating surface state can be formed 
by introducing additional SCOP, such as is or id xy -w&ve 
component [14]. In the present model (d x 2_ y 2 ±is)-wave 
state may occur, because the attractive interaction ex- 
ists between nearest-neighbor sites on a square lattice. 
(It may be possible to induce id xy SCOP near the mag- 
netic impurity with spin-orbit interactions or due to the 
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(11) 



where U = tx if i corresponds to the interface site, U=t 
otherwise. 

The transition between the (d±is)- and (d + p x +p y )- 
wave SC states is of first order, and we did not find a state 
which carries a spontaneous spin current for F /S(d x 2_ y 2) 
interface states. Our results imply that both the state 
with m 7^ 0, A s = and A px (y) 7^ 0, and the state with 
m = 0, A s 7^ and &- px ( y ) — correspond to the local 
minima of the free energy, and that the former (latter) 
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has lower energy when tx > (tx)cr {tx < {tx)cr) with 
(t T )cr ~ 0.02. Actually for t T ~ (i T ) cr , we obtained two 
kinds of solutions depending on the initial values of OP's 
for the iteration, and the energies of two solutions cross 
at tx = (ir)cr- 

V. SUMMARY 

We have studied the proximity effect between uncon- 
ventional superconductors and ferromagnets. It is found 
that a spontaneous spin current can flow along the in- 
terface between a (p x ± ip y )-wave superconductor and 
a ferromagnet due to the coexistence of the magnetiza- 
tion and the chiral superconducting state. In the case 
of a d x 2_ y 2--wave superconductor the induced magnetiza- 
tion destroys the T-breaking (d ± is)-wave surface state 
for relatively large transmission (tx ^ 0.02). For small 
transmission (tx ^ 0.02) the (<i±is)-wave SC state with 
m = is realized, and we did not find the state with 
a spontaneous spin current for any value of tx- This 
implies that only a spontaneous charge current may be 
possible for high-T c cuprates, while both spin and charge 
currents may be expected for spin-triplet superconduc- 
tors, e.g., Sr2Ru04. For conventional s-wave supercon- 
ductors based on the electron-phonon mechanism, the 
ferromagnetism plays a very strong role to suppress the 
superconductivity. On the contrary, in the case of uncon- 
ventional superconductors these effects may be weaker, 
since the superconductivity is caused by the magnetic in- 
teractions. Then we expect it is possible to obtain a sur- 
face state with a spontaneous spin current experimentally 
using spin-triplet superconductors. 
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Fig. 1 Schematic descriptions of the interfaces between 
ferromagnets and superconductors for (a) [100] and (b) 
[110] surfaces. 

Fig. 2 Spatial variations of the magnetization m and the 
SCOP's for a F/S(p x + ip y ) bilayer system with tx = 1, 
U F = 12, U s = 1-5, V s = -2.5, n = -1.6, N x = 60 + 60 
and N y = 120. (a) m, (b) ReA px , (c) lmA py , (d) ReA d 
and (e) ReA s . Note all OP's are non-dimensional, and 
x = corresponds to the interface site of a ferromagnet. 

Fig. 3 Spatial variations of m and the SCOP's for a 
F/S(p x + ip y ) bilayer system with tx = 0.1. Other pa- 
rameters are the same as in Fig. 2. (a) m, (b) KeA px , (c) 
ImApy, (d) Re Ad and (e) ReA s . 

Fig. 4 Spatial variations of m and the SCOP's for a 
F/S(p x + ip y ) bilayer system with Us = 5. Other pa- 
rameters are the same as in Fig. 2. (a) m, (b) ReA px , (c) 
ImApy, (d) Re Ad and (e) ReA s . 
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Fig. 5 The charge and the spin current (in units of 
2eta/h) for F /S{p x +ip y ) bilayer systems, (a) Uf = 12, 
U s = 1.5 and t T = 1, (b) U F = 12, U s = 1.5 and 
t T = 0.1, (c) [Tp = 12, U s = 5 and t T = 1 and (d) 
Uf = 16, f/s = 1.5 and = 1. Other parameters arc 
the same as in Fig. 2. 

Fig. 6 Spatial variations of to, ReA px and ImA py for 
F/S(p x + ip y )/F trilayer systems. The system size is 
N x = 40 + 40 + 40, N y = 120, and other parameters are 
the same as in Fig. 2. (a) to with parallel configuration, 
(b) m with antiparallel configuration in two F layers, re- 
spectively, (c) RcA px and (d) ImA py . Note ReA px and 
ImApy do not depend on the configurations of to. 

Fig. 7 The charge and the spin current (in units of 
2eta/%) for F/S{p x + ip y )/F trilayer systems. Parame- 
ters used are the same as in Fig. 6. Magnetization of two 
F layers are antiparallel in (a), and parallel in (b). 



Fig. 8 Spatial variations of m and the SCOP's for a 
F/S(d x 2_ y 2) bilayer system with tx = 0.1, Uf = 10, 
U s = 1.5, V s = -1.5, n = -0.5, N x > = 60 + 60 and 
N y , = 120. (a) to, (b) ReA d and (c) RcA px (= RcA py ). 

Fig. 9 Spatial variations of to and the SCOP's for a 
F/S(d x 2_ y 2) bilayer system with tx = 0.03. Other pa- 
rameters are the same as in Fig. 8. (a) to, (b) ReA^ and 
(c) ReA px (= ReA P y). 

Fig. 10 Spatial variations of m and the SCOP's for a 
F/S(d x 2_ y 2) bilayer system with = 0.01. Other pa- 
rameters are the same as in Fig. 8. (a) to, (b) RcA^ and 
(c) ImA s . 

Fig. 11 The charge current (in units of 2eta/Ti) for 
F /S(d x 2_ y 2) bilayer systems with (a) tx = 0.1 and 
tr = 0.03, and (b) tx = 0.01. Other parameters are 
the same as in Fig. 8. 
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